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The recently derived nonlocal quantum kinetic equation for dense interacting Fermi systems com- 
bines time derivatives with finite time stepping known from the logistic mapping. This continuous 
delay differential equation equation is a consequence of the microscopic delay time representing the 
dynamics of the deterministic chaotic system. The responsible delay time is explicitly calculated 
and discussed for short range correlations. As a novel feature oscillations in the time evolution of the 
distribution function itself appear and bifurcations up to chaotic behavior occur. The temperature 
and density conditions are presented where such oscillations and bifurcations arise indicating an 
onset of phase transition. 
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The relation between microscopic chaos and ki- 
netic theory within the time evolution of averaged 
distribution functions is a matter of ongoing de- 
bate. If the coarse graining of phase space due to 
averaging is large enough, any microscopic signal 
of chaoticity like bifurcation is lost. The time evo- 
lution of the distribution function evolves conti- 
nously and smoothly in time. If the coarse grain- 
ing is reduced by considering more microscopic 
fluctuations the time evolution of the distribution 
function itself can reveal signs of the underlying 
chaotic motion. To understand this relation dif- 
ferent attempts have been made in the literature 
starting from standard mapping and deriving an 
appropriate kinetic description in terms of finite 
time step equations. Here the opposite route is 
proposed: Starting from a quantum statistical ap- 
proach we employ the nonlocal kinetic theory for 
dense Fermi systems to show that indeed the re- 
sulting kinetic equation can exhibit a nonmono- 
tonic time evolution in a specific temperature and 
density range. This is a consequence of the con- 
sidered many-body correlations. Within a first 
approximation of relaxation times we show that 
a delay differential equation appears which inter- 
polates between finite time stepping like logistic 
mapping and continuous time evolution by differ- 
ential operators. It is argued that the occurring 
bifurcations in the time evolution are a signal of 
onsetting phase transition. 



I. INTRODUCTION 

The kinetic theory of interacting quantum or classical 
gases and the description of deterministic chaotic systems 
are mostly developed separately with very few overlap- 
ping. Almost exclusively for the case of Lorentz gas the 
interlink has been analyzed [1,2], and citation therein. 
The kinetic theory understands the increase of entropy 



as a result of many random collisions. The theory of 
deterministic chaos on the other hand represents the ir- 
reversibility by the characteristic measure of Lyapunov 
exponent. While the kinetic theory can be easily ex- 
tended to quantum systems the quantum chaos is still a 
matter of debate about even the correct definition. 

If both approaches describe some facet of irreversibil- 
ity it should be possible to give relations between them. 
For instance, one can connect the transport coefficients 
with the Lyapunov exponent [3-8]. In [3,6] the fact, that 
the spreading of a small phase space volume is given by 
the sum of Lyapunov exponents, is used to give a rela- 
tion between Lyapunov exponents and viscosity. In [4,5] 
the relation between transport coefficients and Lyapunov 
exponents was presented in terms of Helfand's moments. 
The interlink was possible to establish by reinterpretation 
of the Helfand's moments as stochastic quantities such 
that the mean variance of the time derivatives represents 
just the transport coefficients. In [7] the authors derive a 
density expansion of the largest Lyapunov exponent for 
hard sphere gases from a generalized Lorentz-Boltzmann 
equation. The intimate relation between transport coef- 
ficients and dynamical quantities like the Lyapunov ex- 
ponent is even apparent in quantum Fermi gases in that 
an additional chaotic process behaves like an additional 
relaxation time [8]. 

Recently another more direct relation between kinetic 
theory and deterministic chaotic systems has been estab- 
lished [9] . It appears for the standard mapping that the 
phase-space averaged kinetic description can be formu- 
lated and solved in terms of finite time step equations 
instead of differential ones. We want to address here the 
opposite route starting from a continuous kinetic equa- 
tion for strongly correlated quantum Fermi systems and 
will ask whether there is a reminiscence of the underlying 
chaotic motion in the time evolution of the distribution 
function itself. For a system of hard-sphere gases the ap- 
propriate kinetic description is the Enskog equation con- 
sisting of nonlocal off-sets in the collision integral. Here 
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we will consider the opposite extreme of particles inter- 
acting in such a way that they form short-living molecules 
as correlated states. We will show here that such colli- 
sion duration leads to delay differential equations. In 
particular, it will be found that the time evolution of the 
distribution function itself undergoes oscillations up to 
chaotic behavior similar to Hopf bifurcations observed in 
mappings including memory [10]. This is due to a com- 
petition between attractive potential and repulsive Pauli- 
blocking. It can be considered as a dynamical signal of 
the onset of phase transition due to the formal similarity 
between mean-field phase transitions and bifurcations of 
one-dimensional discrete maps [11]. 

The understanding of the route to chaos has been 
paved with the logistic mapping [12] 



x k+1 = ax k (l - x k ) 



(1) 



which starts to show successive bifurcations for increasing 
parameters 3 < a, the chaos occurring at a = 3.5699... 
An interesting extension of this logistic model has been 
studied by Berezowski [13] in that a finite inertia has 
been added 

+ f{t) = af(t Td)[l f{t Td)] (2) 

with a delay time t<j. For a = this model reduces to 
the logistic mapping (1). In general such retarded dif- 
ferential equations (2) are very complex and are used to 
model periodic and aperiodic dynamics of physiological 
systems [14]. It is possible to simulate nonlinear dynami- 
cal systems by linear delay-differential equations in load- 
ing information in the initial conditions [15]. 

As done in [13] we perform a a linear stability analy- 
sis of (2) according to / = / + Sf with the fixed point 
solution of — /o + a/o — a/o = 0. Without loss of gen- 
erality, we look only at the instable one / = SL ^. The 
Fourier-transformed disturbance Sf(t) = J j^e~ luJt f(u)) 
determines then the complex frequency lu = Q + iT which 
gives for T > the instable and for T < the stable so- 
lutions. One gets from (2) 



-ioja + 1 



ae 



f (l-/o)=0. 



Separating real and imaginary part of (3) one has 

1 + To- = e~ rTd (2 - a) cos n Td 
-Oct = c" rTd (2 - a) sinfh- d . 



(3) 



(4) 



Dividing both sides of (4) gives the equation for the fre- 
quency Q 

Oct + (1 + IV) tan fir d = (5) 
and adding both squared sides of (4) gives 

+ + W<t 2 = c- rT «\2-a\. (6) 



Now we look for instable modes T > which leads to 
y/(l + Ta) 2 + Q 2 a 2 >\2-a\ for r d < 0. (7) 

In the final form we assume Ta <C 1 and we see that the 
presence of the term a shifts the onset of bifurcations and 
of chaos towards higher (lower) values of a if compared 
to the logistic map 



\a-2\ < yjl + (ail) 2 



(8) 



for positive (negative) delay times T d respectively. Fur- 
thermore an oscillation appears with a frequency f2 de- 
termined by (5) 



TdQ + arctan(crO) = n. 



(9) 



There exists a limiting value for the delay time Td below 
which all oscillations dissappear. These oscillations and 
the finite delay time Td are the reasons why Hopf bifurca- 
tions can occur in the model (2) analogously to the model 
considered in [10] where an additional memory has been 
added to (1). 



II. NONLOCAL KINETIC THEORY 

It is now noteworthy to see that the model (2) can in 
fact be derived from microscopic quantum statistics and 
as such bears direct physical relevance. We will show 
that indeed for an interacting Fermi system at low tem- 
peratures the one-particle distribution function obeys a 
kinetic equation which can be written in the form of (2). 
In this way we will express the parameters a and a by 
characteristic physical quantities of the system, the tem- 
perature T, density n, scattering length ao and range of 
the potential r . 

First let us recall the quantum nonlocal kinetic equa- 
tion for the quasiparticle distribution / [16,17] with 
quasiparticle energy s 



dfi dsi dfi 
dt dk dr 



dei dfi 
dr dk 



(l-/l)(l-/ 2 ")/ 3 "/4 



(27T) 5 fi 7 " 
/l/ 2 " (1-/3") (1-/4") 



(10) 



The superscripts — denote the signs of non-local correc- 
tions: h = f(k,r,t), / 2 " = f(p,r-A 2 ,t), / 3 " = f(k-q- 
A K ,r-A 3 ,t-A t ), and f 4 = f(p+q-A K ,r-A 4 ,t-A t ). 
The first (second) part of the collision integral is called 
in(out)-scattering since it (de)populates the distribution 
fi(k,r,t). The time-space symmetries and particle-hole 
symmetries of (10) are discussed in [18]. The form (10) 
used in the present paper, shows explicitly the particle- 
hole symmetry since the out-scattering is the particle- 
hole mirror of the in-scattering. Let us here remark 
only that it can be given an equivalent form with ex- 
plicit time-space symmetry [18]. The scattering measure 
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p = |rr| 2 5(e fc +s p -e k+q -e p _ q -A E )is given by the 
modulus of the scattering T-matrix which is a complex 
quantity T r = \T r \expi<fi. The latter one can be found 
from scattering theory or experimental phase shift anal- 
ysis. All corrections, the A's, describing the non-local 
and non-instant collision [16,17] are given by derivatives 
of the scattering phase shift (f> = Im lnT r (w, k,p,q,t,r) 



£l+£2 



Ae = — 



2 dt 



A 2 - 



A, = 



dp dq dk 



S1+S2 



S1+S2 



dk 



AK = 2-d-r 



S1+S2 



dk 



dcj>_ 
dq 



£l+£2 



(11) 

The nonlocal kinetic equation (10) covers all quantum 
virial corrections on the binary level and conserves den- 
sity, momentum and energy including the corresponding 
two-particle correlated parts [17]. The kinetic equation 
(10) allows for a simple classical interpretation of the 
complicated correlations covered by collisions. The colli- 
sion integral can be understood as a collisional scenario 
that two particles are approaching until they reach their 
correlation distance A 2 then they form a correlated pair 
with the collision duration of A t . During this correlated 
travel they can rotate so that they break up into single 
particles after the time A t at the endpoints A 3j 4. 

The classical kinetic theory of dense gases of Enskog- 
like equations as well as the Landau theory of interacting 
Fermi systems are limiting cases of this nonlocal kinetic 
equation. It requires no more technical problems than 
solving the Boltzmann equation [19,20]. From the for- 
mal point of view this kinetic equation has been derived 
taking into account all terms up to linear order in the 
quasiparticle damping. 

We will now derive the equation (2) from (10) and will 
calculate the required parameters for a model T-matrix 
of separable potential. First we analyze the collision in- 
tegral and assume quadratic dispersion relation for the 
quasiparticle energies. The <5-function in V representing 
the energy conservation in the collision integral reads 

k 2 p 2 (k-q) 2 (p + q) 2 



S< "2m c 
= 5 



+ 



<1\ 



P 

2nib 2m a 2nib ' 

,_k P_, q_ _ <f_,___ + J_x 

y m a m b \q\ 2 m a m b 

S(\q\) 



S(q 0) 



(12) 



with effective masses of particle a, b. The second part 
is the usual considered 5-function excluding zero trans- 
ferred momentum q = 0. This collision part will be 
treated in relaxation time approximation oc (F — f)/r 
where the equilibrium distribution Fa is specified later. 



More interesting is to notice that the q = channel 
leads to an additional part absent in usual local kinetic 
equations like the Boltzmann kinetic equation. To con- 
vince the reader of this novel observation let us rewrite 
the Pauli-blocking factors of (10) for the q — channel 
where f 3 — > /i and f± — > f 2 according to (12) 



(l-/l)(l-/ 2 -)/ 3 "/4- - /l/ 2 -(l-/ 3 ")(l-/4-) 

= / 2 -(i-/ 2 -)[/r-/i 



(13) 



We see that for the case of local kinetic equations without 
delays this specific channel disappears since /{" = f\. 
Therefore we can conclude that in contrast to the usual 
local kinetic equation like the Boltzmann equation, the 
nonlocal equation possesses a hnite zero angle channel in 
the collision integral which is of course of mean-field type 
since no energy or momenta is exchanged. We are now 
going to exploit this observation in a specific way. 

Since the signs of the shifts in the out-scattering parts 
are arbitrary [18] we don't include them at all. Concen- 
trating on the in-scattering part we rewrite 

/ 3 "/4"(l-/ 2 ")(l-/l)l ? =0 

= /r(/r-/i)/r(i-/ 2 -)l fl =o+/3-(i-/r)/r(i-/ 2 -)l <F * 
= /r(/r - /i)/ 2 -(i - / 2 -) + /r(i - /r )/ 2 -(i - / 2 -). 

(14) 

The first part can be expanded 



/f - h 



df a df 



1 dt 

x d f 
-A 3 — 

dr 



' dr 
df 
dk 



A K 



dk 



(15) 



where in the last step we have replaced the time deriva- 
tive of / by the free drift motion = -^x^ + ^x^ 
leading to derivative of on-shell shifts (11). Therefore we 
can absorb the first part of (14) into the left (drift) side 
of the nonlocal kinetic equation (10). These gradients are 
absent for the now considered homogeneous case. Using 
thermal averaging of the occurring microscopic quanti- 
ties r, f and At denoted by < ... > the second part of 
(14) leads to an additional term in the kinetic equation 

df(k,t) f(k,t)-F /(M-A t )[l-/(M-A t )] 



dt 

with 

1 2irs 

r = < — 
xS(e k 

1 2irs 

f =< — 



+ 



(16) 



dpdq 
{2nhf 



\r r (k, P ,q)\' 



£ k+q 



-p-q 



)/ p (l -/*_,)(! -/p+,) > 



/dpdq 
(27Tfi)' 



jiS ^(k,p,q)\ 2 



x6((k-p).±)f p (l-f p )> 
1 d 

A t = -- < — Im In T r (k,p,q, u) \u>=e k +e P > ■ 



(17) 
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The local equilibrium distribution _Fo — /o — t/'?/o(1 — /o) 
is specified here in such a way that equilibrium /o is a so- 
lution of the stationary problem. Conservation laws are 
enforced in the usual way as for relaxation time approxi- 
mation constructing the local equilibrium appropriately. 
We remark that for Fermi ground state /o = (0, 1) we 
have Fo = (0, 1). For the reason of simplicity we will re- 
strict to states above Fermi level F w 0. All discussions 
can be performed with any given F . 

One can see that the equation (2) appears from (16) if 



T 

a = — 

T 



Td 



A*. 



(18) 



Equations (16) and (17) are the main results of the pa- 
per. It shows that from the nonlocal kinetic equation for 
dense interacting Fermi systems one can derive an evolu- 
tion equation which shows the onset of chaos in that an 
oscillation occurs and with increasing parameter a, see 
(8), bifurcations appear and deterministic chaos sets in. 
This is due to the fact that besides the stationary solu- 
tion f of (16) which corresponds to the stationary fixed 
point of (2) we have a branching fixed point which cor- 
responds to 1 — 1/a. For the occupied states F ~ 1 we 
would get the fixed points (1, —1/a) of (2) which would 
result into a minus sign in front of the 2 in equation (8). 



III. MODEL CALCULATION 

To be specific we will now calculate the parameter 
(18) for an interacting Fermi system. We use as an ex- 
ploratory model the separable interaction [21,22] which 
is written in terms of the difference momenta of incoming 
Pi and outgoing particles pi 



V(p!,p 2 ) 



(2irh) 3 \ 



m([3 2 +p 2 )(f3 2 +p 2 ) 



(19) 



in terms of two parameter, the coupling constant A and 
the inverse potential range (3. With the help of this 
potential Bethe-Salpeter equation or retarded T-matrix 
equation T r = V + VG r T r can be solved 



r r (fc,p,o,w) 



(27rft) 3 A/m 



(/3 2 + (k - p) 2 /4)(/3 2 + (k - p + 2q) 2 /4) 

7T 2 A 



x(l 



Calculating the 
cot S = RcT r /ImT r 



f3{Vmhw + i(3) 2 ' 

scattering phase 



(20) 



shift 



pcoi5 = ^)T + ^2r--- a + Yn p +" (21) 

shows that the parameter (3 and A are linked to the scat- 
tering length a and the potential range r in the follow- 
ing way 



ao 



2ft 
J 



(1 



7T 2 A 

4ft 



3ft •*,<, 
r0 = 7 (1 ~W } - 



(22) 



Therefore we can describe with this simple rank-one sep- 
arable potential the scattering length ao and the range r 
of the interaction. Alternatively to the potential range 
one could fit to the bound state energy ui = Eb < which 
yields A = (3(/3 — mE 'b) 2 / 'k 2 . We neglect here medium 
effects on the T-matrix which can be considered as well 
[22]. Now we can calculate the parameter (17) explicitly 
in the low temperature limit. The angular and energy 
integrals can be separated in the usual manner [23,24] 
and we obtain in lowest order temperature T 



7rT 2 



a 2 s ci(P/p f )i 



c 2 {(3/p f )s 



t 3fte/ 
1 4T 2 (3 2 

T 7rft pj 

A * = -^c 3 (a s ,P/p f ) 



(23) 



with the Fermi energy e/ = p"j/2m and the dimension- 
less scattering length a s = aop / /ft. The functions c\ , c 2 
describing finite potential range effects and read 



(x 2 + 1 - y 2 ) 4 



2y/ 2x e 1 

2 ) 3 V^y ~ x 2 



c 2 (x) = [ dy- 

J {l-y+2x 
-l 

ax 4 (-6 + 3aa; + 2a; 2 



l-- + o(x- 4 ) 



c 3 (a,x) = 



(l + x 2 ) (4.x 4 + (2-ax) (2-ax + 4:X 2 -4:ax 3 )) 



4(1 + a 2 ) 



+ o(x 



(24) 



where the first two integrals are elementary functions. 
The final result for the parameter of (2) are now rep- 
resented in terms of the microscopic scattering length, 
a s = pfa /h, and the range of the potential (22) as 

_ e f 12c 2 (3 2 
a ~ T ir 2 Cl p 2 
_ fte/ 3 1 



rp2 



itcis at 



(25) 



From this we can now give the conditions for the onset 
of bifurcation and route to chaos according to (8). 

First it is interesting to discuss the sign of required 
for (8). Since c 3 in (24) determine Td via (23) we obtain 



4 



ao < : T d ^ for p f ^ p c 



0<a Q < f: 
2h 



for pf ^ 



Pc 



(26) 



a > f : T d > 



with p c = P/ \J"i — ^§7p- It shows nicely the interplay 
between the range of potential and the density. For lower 
densities than the critical one, pf < p c , the delay time 
follows the sign of scattering length. For higher densities 
the delay time flips the sign if the scattering length is 
lower than 2h/f3. 

Now we analyze (8) and (9) for the low temperature 
expansion compared to the Fermi energy. The final re- 
sulting condition for the occurrence of oscillations with 
the frequency f2 reads 



T d < 0: T < e f 



3^ 



O — e f ( 31T I g£l ( Pf a ° \ 2 (Z_X 

" — ft I 2c 3 18 V h ) \e f ) 



T d > 0: T > Ef 



(jLY 



(27) 



where C 



e^TSf < 1 accordin g to tne expan- 

sion with respect to the range of the potential. These 
conditions for the temperature together with (26) allow 
to calculate the density and temperature range for which 
the oscillations appear and chaotic behavior is possible. 

In the end we should discuss the validity range of the 
kinetic equation. The nonlocal kinetic equation (10) we 
have started from has been derived under the assumption 
that the delay time is smaller than the relaxation time. 
This translates into r d < a and from (25) we obtain 



T<6f 



12 1 

7TSCiC3 a s 



(28) 



which has an overlapping region with (27) and is con- 
sistent with the low temperature expansion used here. 
To illustrate this we want to consider nuclear matter 
with realistic parameter. The triplet channel has a scat- 
tering length of a = 5.379fm and the deuteron bound 
state energy is reproduced if the parameter for the sep- 
arable potential is = 1.4488/i/fm [21]. Then the 
third case of (26) applies. The corresponding third case 
of (27) together with (28) gives then the criterion on 
the temperature 0.118e/ > T > .0813/le / / 'imp f which 
makes only sense if the density or the Fermi momentum 
Pf > 2.24ft/fm which corresponds to about four times nu- 
clear matter density no = 0.14fm -3 . The corresponding 



temperatures will be around lOMeV and the oscillation 
period 1/51 w 2fm/c. These parameter are reached in 
experimental heavy ion reactions. Therefore we will find 
a parameter range where bifurcation and the route to 
chaos can occur in the time evolution of the distribution 
function itself. 



IV. CONCLUSION 

Summarizing, dense interacting Fermi system are de- 
scribed by nonlocal kinetic equations which can be shown 
to be of the type of delay differential equations. They in- 
terpolate between finite time stepping of the logistic map 
and continuous time derivatives. The important novel 
observation is that the distribution function itself shows 
temporal oscillations around the equilibrium value with 
a basic frequency of about the Fermi energy. For certain 
density and temperature conditions, (27), this oscillation 
bifurcates and the route to chaotic time evolution sets 
in. It could be considered as a signal of onsetting phase 
transition in accordance with the analogy found in [11] 
between mean-field phase transitions and bifurcations of 
discrete maps. 

In the end it should be commented on the entropy. It 
is clear that from the above described possibility of oscil- 
lations in the time evolution of the distribution function 
follows that also the entropy will oscillate. This has been 
seen already in the collisionless but selfconsistcnt Vlasov 
equation [25,26] which appears if we neglect collisions in 
(2) at all and calculated the quasiparticle energy on the 
drift side from the selfconsistcnt mean field. The inclu- 
sion of the collisions are expected to damp this oscilla- 
tions. As we have shown above there might be a regime 
where the correlations are such that oscillations pertain. 
In such cases the entropy would oscillate. The complete 
balance equations [17] which follows from the nonlocal 
kinetic theory (2) show that the observables contain be- 
sides the Landau one-particle parts also correlated two- 
particle parts. Though the explicit forms are known for 
the density, current, energy and stress tensor and their 
corresponding conservation laws are proven [17], we have 
not been able yet to reconstruct the explicit form of the 
entropy beyond the single particle expression. Therefore 
we could not find an explicit H-theorem with two-particle 
correlations included. This will remain to future work. 
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